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Abstract 

We demonstrate the actuation of a double beam opto-mechanical cavity with a sinusoidally varying 
optical input power. We observe the driven mechanical motion with only 200 nW coupled to the optical 
cavity mode. We also investigate the pump power dependence of the radio-frequency response for both 
the driving power and the probe power. Finally, we investigate the dependence of the amplitude of the 
mechanical motion on mechanical cavity quality factor. 

Optoniechanics, the study of the interaction between light and mechanical motion, has recently captured 
the imagination of photonics researchers For example, researchers have probed radio frequency (RF) 

mechanical motion of nanometer sized objects [31 0]. In addition, proposals for using the optical gradient 
force to induce mechanical motion |6l [3 [5] with [H [3 IH [TOl [11] and without [T2j [131 HI] the use of an optical 
cavity have been experimentally demonstrated. In fact, at very high optical and mechanical confinement 
[Tni[Sl[in], the amplitude of a mechanical mode can be greatly increased. In addition, at high input powers, 
regenerative mechanical oscillations occur, where the linewidth of the mechanical mode greatly decreases, 
while the amplitude of the mechanical oscillation greatly increases. 

The experiments above have been done with continuous-wave (CW) excitation of an optical cavity or 
modulated excitation of a waveguide. However, the CW excitation mechanism requires the mechanical 
motion to induce an out-of-phase modulation of the laser input, as only those forces in quadrature with 
the mechanical motion perform mechanical work on the structure. Such effects are generally small, as the 
thermal motion of the structure only weakly perturbs the optical transmission properties of a waveguide or 
cavity. In the CW case, the amount of work done on the mechanical cavity is proportional to where 
K is the optical field decay rate in the cavity. However, an alternative to increase the transduction between 
optical power and mechanical motion is to use modulated pumping |12[ [T5] . Such a scheme can do work 
that is proportional to k~^, and greatly reduce the amount of power needed to excite the mechanical mode 
(see Appendix). In this work, we demonstrate the use of optical pump modulation in conjunction with 
an optical cavity to reduce the amount of power needed to actuate the mechanical mode. Because of the 
optical confinement and recirculation of photons, we hope to obtain large mechanical oscillations without 
regenerative feedback. 

In particular, we choose to work with the double beam one-dimensional photonic crystal (PC) cavity 
configuration in silicon. Due to its high optical quality factor {Q > 10*), which enhances the circulating 
optical power inside the cavity, and low mode volume (~ (X/n)^), which also enhances the local field potential, 
the PC cavity can greatly enhance the optical gradient force. The optomechanical coupling rate is defined 
as: 

duj 

90M = ^^ (1) 

where uj is the optical cavity frequency and x is the mechanical displacement of the cavity. By using cavities 
where the £'-field is increased near material boundaries (such as in a slotted design (HUT]), the frequency 
perturbation with mechanical motion and the optomechanical coupling can both be tailored. 

We fabricate devices on a silicon-on-insulator (SOI) wafer with a 150 nm thick layer of Si and a 1 /im 
thick oxide layer, such as the cavity shown in Fig. [ija). The beam cavities have lengths of approximately 13 



/^m, single beam widths of 550 nm, and a middle slot width of 100 nm. We use the design of Ref. [TB], where 
the hole lattice constant is kept constant at a = 400 nm, and the radii of circular holes are reduced as the 
distance from the center of the cavity increases. The larger holes at the center of the cavity create an optical 
potential well that lies in the optical bandgap of the array of outer holes, and such a design allows robust and 
high-efficiency coupling to the cavity region via a coupling waveguide. The hole at the center of the cavity 
has radius r = 0.28a, and the total cavity length is 34 holes. The cavity is fabricated with electron-beam 
lithography, and the pattern is transferred into the silicon layer by a Cl2:HBr plasma dry etch. The oxide 
sacrificial layer is then etched away using a buffered oxide etch (BOE) to obtain the free standing beams. 
In addition to the beam cavity, we also attach coupling waveguides on both sides of the cavity, and one of 
the waveguides is bent 90° to configure the device to be probed in a cross-polarization geometry [Fig. [Ija)] 

[MUD]. 





Figure 1: (a) Scanning electron microscope image of the fabricated cavity. The Ey field of the (h) TE\^+ and 
(c) TE2,+ optical modes, (d) The first order common in-plane mechanical mode, and (e) the first order differential 
in-plane mechanical modes are plotted with the color map assigned to the in-plane (y ) motion. 
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We first simulate the beam cavities in the optical regime using the three dimensional finite-difference time- 
domain (3D-FDTD) method. Double beam cavities support bonded and anti-bonded (— ) optical super- 
modes, formed from the the transverse-electric (TE) modes of the individual beam cavities. In particular, 
the Ey field is symmetric or anti-symmetric about the xz-plane going through the slot for the bonded and 
anti-bonded modes, respectively. We find that the first (TEi.+) and second (TE2,+) order bonded modes 
(see Fig [Hb)-(c)) have theoretical radiation-limited Qs of 30,000 and 1,500, respectively. We observe an 
enhanced electric field in the air slot region for the bonded modes because of the continuity conditions for 
the dominant Ey field at the slot boundaries (i.e. continuity of the displacement vector eE). Thus, we expect 
that the bonded optical modes have the highest optomechanical coupling to the in-plane mechanical modes, 
as the high electric field concentration in the middle of the cavity enhances the change in the optical cavity 
frequency with mechanical deformations. For this reason, we work with the first and second order bonded 
optical modes in our experiments. 

We experimentally analyze the optical properties of the cavity using the setup in Fig. [SJa). We pump the 
cavities with a broadband LED bank, which is coupled into a waveguide using a dielectric grating coupler. 
We align the cavity such that the input grating polarization is along \H), while the output polarization is 
along \V), to obtain the maximum signal to noise ratio. The transmission characteristics of the cavity are 
shown in Fig. ^h), where we are able to observe the first two orders of the bonded and the anti-bonded 
modes. We are able to differentiate the bonded modes from the anti-bonded modes by moving the input 
beam on the grating coupler to change the input parity. The first order modes have high Q-factors, and 
we use a tunable laser to fully characterize the cavity. The laser scan at low input powers (1 nW) shows a 
Lorentzian spectrum with Q w 15, 000 for the bonded first order mode (TEi_+) (inset Fig. [Hb)). In addition, 
we observe that the higher order bonded mode (TE2,+) has Q k, 2, 000. Both Q values are comparable to 
the FDTD simulated values. 

We next use the COMSOL finite element solver to find the frequencies of the mechanical modes, using 
library parameters for sihcon: Young's modulus of 131 GPa, Poisson's ratio of 0.27, and a density of 2.33 
g/cm^ . As described above and in previous work [4], mechanical modes with in-plane (in this case, referring 
to the xy plane) motion will have significant optomechanical coupling to the bonded modes. In particular, 
we find the first order common and differential modes for in-plane motion j4]. The common and differential 
modes have the beams moving in phase and out of phase, primarily in the y-direction, and have displacement 
profiles shown in Fig. [TJd) and (e), respectively. By simulating the structure observed in the SEM image, 
we find that these two mechanical modes have mechanical frequencies of 25.72 MHz and 26.74 MHz. We 
find the optomechanical coupling strength similarly to previous work [J" , with the optomechanical coupling 
length defined as: 

1 _l !dA[^^-h){^e\E^^\^-A{e-^)\D^\^) 
LOM 2 SdVe\E\'^ ' 

Here, q is the mechanical displacement, a is the parameterized displacement of the mechanical mode, fi is 
the surface normal vector, is the electric field parallel to the surface, D± is displacement field normal 
to the surface, Ae = ei — 62, and A(e)~^ = ej""^ — e^"^, with ei being the dielectric constant of silicon, and 
£2 the dielectric constant of the surrounding medium. Because of the high iJ-field enhancement in the slot 
and the differential mechanical resonance having opposite parity to the Ey field, we observe very strong 
optomechanical coupling lengths of ^qM — 1-3 A'ni ^^nd 1.8 fim for the coupling between the differential 
mechanical mode and the TEi^-|- and TE2,+ optical modes, respectively. On the other hand, the coupling 
between the TEi^+ and TE2,+ optical modes and the common mechanical mode was calculated to be far 
weaker (^QM because this mechanical mode has the same parity as the optical field. 

In order to first characterize the mechanical modes of the system, we pump the second order bonded 
mode with a red detuned probe laser, at the cavity half-max, with low pump power (300 /iW before the 
objective) to observe the mechanical modes in air. The transmission signal is fiber-coupled and sent to a 
photodiode detector with a transimpedance gain of 2.5 X 10"* V/A and a bandwidth of 125 MHz, and the 
electrical signal is then read by an RF spectrum analyzer. We estimate coupling efficiencies of 2% to the 
TE2,+ mode and 0.5% to the TEi__|_ mode, assuming symmetric losses at the input and output gratings, 
and accounting for the transmission losses of the coupling waveguides using FDTD simulations. We observe 
the two mechanical modes in the RF spectrum, shown in Fig. El^a), which correspond well to the simulated 
in-plane mechanical mode frequencies, and slight discrepancies can be attributed to minor differences in 
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Figure 2: (a) The optical setup used to probe the optomechanical cavity, (h) Spectrum of the cavity observed in 
transmission using a broadband LED. The first and second order bonded (-\-) and anti-bonded (—) modes are labeled. 
The inset shows a laser scan of the TEi^+ cavity mode for excitation, with a fit to a Lorentzian Imeshape having 
Q ^ 15, 000. 



the clamping conditions of the fabricated device. Because of the low optical Q of the TE2,+ mode and the 
low optical power buildup, we do not observe the giant optical spring effect seen in previous works [4l ITT], 
as the mechanical modes do not change frequency with increasing pump power. We also do not observe 
significant changes in the optical cavity wavelength with pump power, suggesting minimal heating. Because 
of the mechanical damping of the ambient atmosphere, the mechanical Q-factors of these modes are limited 
to 50-100. When we test the same cavity in vacuum, we observe the two modes more clearly, as shown in 
Fig. [3ljb). In vacuum, the mechanical Qs are as high as 2,500, and are limited by the clamping geometry of 
our cavity. We choose to work with the higher frequency mode (the differential mode) , as it is the in-plane 
mechanical mode with higher optomechanical coupling to the second order optical mode. 

Next, we pump the TEi__|. mode with a second (pump) laser tuned to the optical cavity resonance 
wavelength and sinusoidally modulated near the RF frequency of the mechanical mode, while keeping the 
first CW laser tuned to the half-maximum of TE2,+ . We observe the effect of the second, modulated laser on 
the RF modulation of the first laser. We scan through the first-order optical mode with various unmodulated 
powers, and observe that the first order cavity resonance is not significantly changed, suggesting that the 
injected power on the first order optical mode does not change the temperature of the beam, and thus does not 
modulate the beam transmission via the thermo-optic effect. Although both lasers pass through the cavity 
and are extracted with the same output grating coupler, the laser on TEi__|_ is blocked by a band-pass filter 
centered at 1550 nm with a full-width at half-max of 12 nm. The power of the laser on TEi._|_ is modulated 
by a Mach-Zender interferometer modulator with a bandwidth of 2.5 GHz and full modulation depth [Fig. 
[5]Ja)]. First, we fix the input power on the first-order optical mode at under 2 ^W, and scan the modulation 
frequency through the mechanical resonance. When we tune the RF input frequency near the mechanical 
resonance frequency, we observe a narrow response in the RF spectrum (of the laser on TE2,+) [Fig. EJc)]. 
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Figure 3: The RF spectrum of the mechanical modes under study in (a) ambient atmosphere, and in (b) vacuum, 
(c) The time averaged spectrum of the differential mechanical mode from part (b) is shown (green points), observed 
as RF sidebands of the laser tuned to TE2,+ . The non-averaged RF spectrum showing the sharp RF response when a 
modulated laser on TEi^+ is added is also plotted (blue line). The inset shows the same data zoomed in, to observe 
the thermal driven mechanical mode in the background, (d) The integrated power within the sharp RF response of 
the laser on TE2,+ [from (c)J with different RF modulation frequencies of the laser on TEi^+. The two dotted curves 
correspond to two different average input powers on the first order mode and fixed input power on the second order 
mode. A closer zoom of the mode shown in part (b) of the figure is shown as a reference at the bottom (blue). 



In addition, as the RF input frequency is tuned around the mechanical resonance frequency, we observe 
that the integrated power within the narrow bandwidth response matches exactly that of the mechanical 
cavity resonance [Fig. Eld)], suggesting that the optical power in the first order mode is modulating the 
transmission properties of the second order mode through the mechanical resonance. In addition, we observe 
the Lorentzian mechanical mode with far better signal to noise, and can observe the tails of the mechanical 
mode even when detuned by more than three mechanical cavity linewidths. 

We also measure the RF response of the probe laser on TE2,+ as we change the power of the modulated 
pump laser on TEi,+. We first do so with the probe power for TE2,+ fixed at 2 fiW coupled into the cavity, 
and observe the RF response with varying average power on TEi,+ for different RF detunings from the 
mechanical resonance [plotted on a log- log scale in FigHl^a)]. Similar to the data in Fig. [3jd), we observe the 
RF response is decreased as the modulation frequency is detuned from the mechanical resonance. We observe 
that the relationship between the integrated power in the RF response and the input laser power on TEi.-|_ is 
quadratic for all detunings. This is expected, as the RF spectrum analyzer measures the power of the voltage 
signal from the transimpedance amplifier of our detector, and that power has a quadratic relationship with 
the amplifier output voltage and thus a quadratic relationship with the output RF oscillation amplitude. 
This indicates a linear relationship between displacement and input pump power on the first order mode. 

We also measure the RF power spectrum from TE2^+ when we fix the average laser power on TEi__|_, 
and increase the power of the pump on the second order mode, as shown in Fig|4l[b). Again, we observe 
that the integrated RF response of the driven mechanical mode is quadratic with the input power, which is 
expected as the sideband amplitude is linearly related to the probe power. We also obtain the RF response 
as a function of the input power on TEi^+ for various probe powers on TE2,+ , shown in FigUJc). The 
RF response is reduced for lower input powers, as the sideband powers are proportional to the input probe 
power. However, we are able to observe an RF response with only 100 nW coupled to the TEi^+ mode 
to drive the mechanical oscillations, and only 200 nW coupled to the TE2,+ mode to sense the mechanical 
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Figure 4: (a) The integrated intensity in the RF response collected from TE2,+ as a function of average input power 
on TEi,-f- for different detunings of the RF modulation frequency from the mechanical resonance at a fixed probe power 
(2 i-iW) on TE2.+ . (b) The integrated intensity in the RF response as a function of different probe powers on TE2, + , 
at two different fixed average pump powers on TEi^+. (c) The RF response as a function of input power on TEi^ + 
with different probe pump powers on TE2.+ . (d) The integrated RF response as a function of average pump power on 
TEi^^. The two curves correspond to the response at ambient atmosphere and in vacuum, both with the same probe 
intensity on the TE2,+ mode (2 nW). 

motion. 

Finally, we compare the efficiency of exciting the mechanical mode in vacuum and in ambient atmosphere. 
We fix the input power for the probe laser on the TE2.+ mode in both air and vacuum to 2 /iW, and obtain 
the same output coupled power into our photodetector. We obtain the power series from the same cavity 
under both conditions, which is shown in Fig. lUJd). As expected, the amplitude of the mechanical oscillation 
is significantly higher in vacuum than in ambient atmosphere, due to the higher mechanical Q. In fact, the 
experimentally measured factor of 20 between the power needed to generate the same RF response in air 
and vacuum matches well with the ratio of mechanical Qs for the two conditions (31). 

In conclusion, we have demonstrated resonant actuation of a mechanical mode with optical gradient 
forces. The input power needed to observe driven motion of the mechanical cavity is greatly decreased in the 
presence of an optical cavity, and hundreds of nanowatts can drive the mechanical motion via a modulated 
laser coupled to a second cavity mode. This type of excitation can be used to probe various mechanical modes, 
as the RF response can be increased relative to the thermal-driven oscillations. Furthermore, optomechanical 
cavities can be used to mix RF signals, with the mechanical resonance enhancing the beat note of two RF 
signals. Similarly, the actuation of mechanical motion can also be used for a variety of applications, such as 
mechanical motors that do work on nanometer-sized objects. 

We acknowledge support from the Presidential Early Career Award for Science and Engineering (PECASE), 
administered by the Office of Naval Research (Dr. Chagaan Baatar). We also acknowledge support from 
the National Science Foundation graduate research fellowship (YG), and the Stanford Graduate Fellowship 
(AR, AM). 
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1 Appendix: Theory of resonant excitation of mechanical mode 
with optical gradient force: 

We would like to solve for the mechanical amplitude as a function of the average input power of a modulated 
laser. We follow the derivation given in Ref. [3 and start with the cavity field equation: 

ia{t)ujQ ^ [k, 



c(t) = - (f + ^.o) c{t) + '-^c{t) + J^sit)e-' (3) 



where s{t) is the time-varying pump field, luq is the cavity frequency, k is the cavity field decay rate, Ke 
is the external coupling rate, Lqm is the optomechanical coupling, c(t) is the cavity field, and a{t) is the 
mechanical mode amplitude. In this case, we are inputing a laser at uj which is detuned from the optical 
cavity mode center frequency, and the input is modulated periodically with frequency fl, which is detuned 
from the mechanical mode center frequency fio- 

We assume sinusoidal mechanical motion, such that the beam also moves with modulation frequency fl: 

a[t) = ao sin(m) (4) 

Note that Vt could be different from VLq, but since we're driving the motion, we can assume the mechanical 
mode responds with the same frequency. Then the equation becomes: 



2 



The homogeneous solution is: 



Ch{t) = Co/u = Co exp + iwojt 



2 LoM^ 
= Co exp (-(| + iu;o)t) ^ (-z)" J„(/3)e"'"* (6) 

n 

with /3 = ujQao/ Lqm^, and the inhomogeneous solution is: 

Since our pump is modulated with frequency f2, we express s{t) = ^ afce**"'^* as a Fourier Series, and find 

k 

the full inhomogeneous solution: 



n ' fc 



_ [ — iA.-\-i[n-\-k)Q,~iujQ)t — i/3 cos(£7i) /q\ 

q-iA + i{n + k)n ^' 



.,k 2 



with A = ui — ojQ, and neglecting the Ke/2 term as normalization: 

Because the homogeneous solution levels out with rate k and this is fast, the particular solution is the 
steady state solution. The optical force is: 

kp(^)P ^ 1 i"'~™Jnil3)Jm{l3)akai ^i[{n+k)-{m+l)]nt rQ\ 

LoM LoM A'^-iA + i{n + k)VL){!^+iA-i{m + l)n) ^' 

Taking only the zeroth order in Jo(/3), as /3 ^ 1 and Ji(/3) w (3: 

\cp{t)\^ _ J§{(3)akat t{k-i)nt .-.(^n 

k,l 



LoM ioM ^ (f - iA + ifcfi)(f +iA - iZf2) 
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1.1 Example 1: Cosine input 

Let's input s{t) — so{l + cos{flt))/2: 



sit) = so U + 



(11) 



Thus we have aq = 1/2, ai = a_i = 1/4. 

The normahzation for the time dependent portion of the input is J |s(t)P — J{1 + cos(ili))^/4 = 
{Sir) /(An). We want to keep the average power the same, so A'^{3n)/{4n)/T = 1 = A2(37r)/(4ri)/(27r/17), 
or A = ys/S- 

The optical force, normalized to the average input power is then: 



IsopKe^^ LOM ^ ( 



A + ikn){i + iA - iin) 



(12) 



We wiU only consider the elements with frequency fi, in quadrature with the beam motion, as they will 
contribute to work getting done on the mechanical mode, so we isolate the cos(f2i) terms: 



F 



4L 



OM 



cos(rii) 



1 



f + (A - f7)2 f 



AO. 



[A + ny 
An 



(^ + A2)(f + (A - nr) + A^){^ + (A + nf) 



(13) 



Note that this force is maximized near A = (as all four terms are near Lorentzian functions in terms 
of A), and we consider the force amplitude (dropping the harmonic variation): 



F 



sol 



jm 

'■KeA^ 2LoM 



2LoM 



or with the normalization (such that input power is proportional to [sp): 

F 8 1 1 



We note that in the case of pumping with a CW laser, the equivalent force is [3]: 



(14) 



(15) 



F 
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IsoP^e 2LoM 



4KAn^ 



(f + A2)(f + (A - mi^ + (A + m 



(16) 



Thus, comparing some sort of AC pump scheme (assuming A = 0, to maximize force) to the DC pumping 
(assuming A = k/2, where the force is approximately maximized), we see that the transferred power should 
be approximately K^/(/3f2^) more efficient. In addition, if we assume that our in-coupling efficiency is 
sufficiently high, then we would have ^ k. Using the above two equations, our optical force is oc k^^ for 
the CW case, and oc for the modulated laser case. 

Note that our thermal amplitude is (x^) — kbT/m^^fl^, and « 10 pm in this case, which places us 
in the high /? regime (despite the sidebands being unresolved, we have (3 — 45). We can calculate the force 
as a function of (3 as well, plotted in Fig. \5\ using real parameters of nieff =2x 10~^^ kg, i7o = 27r-22 x 10^ 
Hz, Qm — 70, Ke = k/2, Lqm —2 fim, optical wavelength A = 1500 nm, and optical Q ^ 2 x 10**. We 
plot the kernel of the force term using Eqn. [S]for different /3, but fixing all other parameters, and plot the 
results in Fig. [5] We notice that for our parameters, the force on the beam is relatively unchanged even up 
to /3 w 100. Thus, we use sinusoidal pump to increase the force amplitude. 
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Figure 5: The theoretical average force on the mechanical mode for a fixed average input power of the modulated 
input, as a function of p. 



1.2 Example 2: Square wave input 

We can also explore input powers that are periodic with frequency f2, but not sinuisoidal. One example is a 
square wave input that is that is zero for some amount of time, and a fixed amplitude A = ^ (chosen to 
have fixed energy input) for time 2Ti. The pulse train is input with period T = 27r/f2. 

^^^^ 1 otherwise ^^'^^ 




The Fourier coefficients of this input are: 



j 2ATi/T if A: = 

^sin(.^T.) ^^^^^^.^^ (18) 



kn 
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By evaluating the sum numerically, we obtain the force as a function of the duty cycle {2Ti/T) in Fig. 
[6] for A = (other detunings only decreased the force) . We observe that the forcing term that does work is 
not drastically increased with pulsed (short duty cycle) pumping. 

Note that this is assuming that only the zeroth order correction for J„(/?) is necessary. It is possible 
that higher order corrections, like that used in the derivation from [3], may be needed. We observe that 
the maximum average force resulting from a square wave input is lower compared with the sinusoidal input 
of the same average power, which is expected as the power of the square input is spread into more Fourier 
components. 
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